Configuration space in electron glasses 






CD 



X3 

o 



> 

o 
o 

ON 

On 
-(— > 



o 



X 



A. Perez-Garrido, M. Ortuiio, A. M. Somoza and A. Diaz-Sanchez 
Departamento de Fisica, Universidad de Murcia, Murcia 30.071, Spain 

We study numerically the configuration space at low energy of electron glasses. We consider 
systems with Coulomb interactions, short-range interactions and no interactions. First, we calculate 
the integrated density of configurations as a function of energy. At a given energy, this density is 
smaller for Coulomb glasses than for short-range systems, which in turn is smaller than for non- 
interacting systems. We analyze how the site occupancy varies with the number of configurations. 
Through this study we estimate the number of particles involved in a typical low-energy transition 
between configurations. This number increases with system size for long range interactions, while 
it is basically constant for a short-range interaction. Finally we calculate the density of metastable 
configurations, i. e. valleys, classified according to their degree of stability. 



PACS numbers; 71.10-w, 73.61. Jc 
I. INTRODUCTION 

The term Coulomb glass refers to Anderson insulators 
with Coulomb interactions between the localized elec- 
trons, and it has established itself as an important model 
in the study of the electronics properties of insulators [Q . 
It assumes that quantum energies t arising from tunnel- 
ing are much smaller than the other important energies 
in the problem, i.e. Coulomb interactions and random 
energy fluctuations. In this situation the Hubbard energy 
is very large and practically forbids double occupation of 
sites. The interesting case then corresponds to a number 
of electrons roughly half the number of sites. 

Until a few years ago, most studies of the Coulomb 
glass consisted mainly of finding the ground state of the 
system, and considering from there on only one-particle 
transitions H-H]. The main problem arising with this is 
the ncglcction of many-body effects. Several experimen- 
tal [||-0 and numerical |^ results give evidences that elec- 
tronic correlations are very important in these systems 
and they must be taken into account. More recently, 
methods were developed to obtain an almost complete 
set of low-lying states of the Coulomb glass [^,0 . This 
process allows for a much more detailed consideration of 
many-body effects. 

A proper understanding of the low-energy configura- 
tion space is extremely important for the study of most 
physical properties. Previous works on Coulomb glasses 
have never focus on the properties and relations between 
configurations. We only have an indirect knowledge of 
them through the analysis of different physical prob- 
lems. In this work we perform a systematic study of the 
low-energy configurations of three-dimensional Coulomb 
glasses. For the sake of comparison, we have also com- 
puted the same properties for non-interacting systems 
and for disordered systems with short range interactions. 

In the next section, we introduce the model used and 
details of the numerical simulations. In section III, we 
study the entropy of the systems and compare to that ob- 
tained by two different analytical models. In section IV, 



we obtain the average difference in occupancy between 
the ground state and configurations as a function of en- 
ergy, and in section V we calculate the average number 
of electrons participating in a typical low-energy tran- 
sition. In section VI we study the density of valleys in 
configuration space. Finally, in section VII we extract 
some conclusions. 



II. NUMERICAL MODEL 

We consider three-dimensional systems deep in the in- 
sulating regime, where the relevant energies in the prob- 
lem are the Coulomb interaction and the disorder energy. 
This corresponds to the condition a <^rQ, where a is the 
localization radius, tq = (47rn/3)~^/'^ is the typical dis- 
tance between sites and n the concentration of sites. 

We use the standard tight-binding Hamiltonian: 



i^ = ^e.K-if)+^F,„ 



(1) 
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Ci are the site random energies, uniformly distributed at 
random in the interval {—W/2,W/2), W being the dis- 
order strength. We take the number of electrons to be 
half the number of sites and K — 0.5. The sites are 
arranged at random, but with a minimum separation 
between them, which we choose to be 0.5ro. We take 
e^/ro as our unit of energy and rg as our unit of dis- 
tance. Vij is the interaction potential. In this work we 
use three types of potentials, coulombic, short range and 
non-interacting. The Coulomb potential is given by the 
expression 
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where r^- is the separation between sites i and j. For the 
short range potential we choose 



V,, 



[n, - K) (n, - K) 



(3) 



G was taken equal to 0.7. Finally, the non-interacting 
case simply corresponds to Vy = 0. 

A numerical algorithm developed to obtain the ground 
state and the lowest energy many-particle configurations 
of these systems has been used for the calculation of low 
energy states pi| . This algorithm considers simultaneous 
many-electron transitions constructed from (of the order 
of ten) independent one-electron transitions and perform 
the combination that relaxes more energy. Starting from 
a random initial configuration, we repeat the relaxation 
procedure until we cannot reduce the energy any more. 
The whole scheme is repeated for different initial ran- 
dom configurations of the charges until the configuration 
of lowest energy is found ten times. The configurations 
thus generated were memorized in terms of site occu- 
pation numbers (0 or 1) and of energy, whenever this 
was less than the highest energy configuration in mem- 
ory storage. We complete the set of low-energy configu- 
rations by generating all the states that differ by one- or 
two-electron transitions from any configuration stored. 



III. NUMBER OF CONFIGURATIONS 

The first magnitude to study is the number of con- 
figurations as a function of energy. Let us call density 
of states to the number of one-particle states per unit 
energy, and density of configurations n(Ei) to the num- 
ber of many-particle configurations per unit energy. In 
particular, we consider the integrated density of config- 
urations N(E), i.e. the number of configurations with 
energy smaller than or equal to E. If the occupation 
probability of all configurations up to an energy E is the 
same, the entropy, S, is then equal to fee In A^(_E), where 
fee is the Boltzmann constant which we take equal to 1. 
The number of configurations strongly depends on the 
size of the sample, but this dependence can be roughly 
taken into account by an adequate normalization of the 
energy. The idea is to measure the energy in terms of 
the average energy spacing between single particle states 
near the Fermi level (in the absence of interactions): 



E = SV/W, 



(4) 



where E is the normalized energy, V the volume of the 
system, W the disorder strength and £ the energy to be 
normalized. This procedure provides a density of config- 
urations independent of dimensionality and system size 
for non interacting systems. 
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FIG. 1. \nN{E) versus energy for systems with no interac- 
tions (solid circles), short range interactions (triangles) and 
Coulomb interactions (empty circles). The solid line corre- 
sponds to the standard result, Eq. (H), and the dashed line to 
our discrete model, Eq. (h_5|). 

In Fig. 1 we plot on a double logarithmic scale the 
average of the entropy as a function of the normalized 
energy, i.e. measured in terms of the one-particle en- 
ergy spacing. Averages are taken over twenty samples. 
The solid dots correspond to the non-interacting case, the 
triangles to a system with short-range interactions, and 
the empty dots to a Coulomb glass. As we have previ- 
ously mentioned, the results for the non-interacting case 
are exactly independent of the sample size. The two in- 
teracting cases correspond to averages over samples with 
a number of sites equal to N^ — 900. Although, in prin- 
ciple, the results for interacting systems depend on size, 
the variation is very small and we have preferred not to 
show them for several sizes for the sake of clarity. In all 
cases considered, the fluctuations from sample to sample 
are large, even in the logarithmic scale considered, but 
averages are relatively well behaved. 

For non-interacting systems it is possible to make the- 
oretical calculations for the entropy. The classical result 
obtained from statistical mechanics for an infinite sys- 
tem with a constant density of states corresponds to the 
solid line in Fig. 1. The entropy is proportional to \/E, 
and so this is a straight line with slope equal to 0.5. 
Although the previous result is well known, we deduce 
it below in order to express it in our notation and be- 
cause the method of calculation can be extended to other 
non-standard magnitudes as we will see later on. The 
dashed line connects the points obtained from a finite 
discrete model that we also describe in the next subsec- 
tion. All systems have a integrated density of the form 
N{E) (X exp{E / Eq)"' . In all cases, for the sizes consid- 
ered we checked that a ~ 0.7. The value of the exponent, 
a, is deeply affected by finite size effect. For interacting 
systems we are not sure of the value the exponent in the 
thermodynamic limit. The similarity between the inter- 
acting and non-interacting systems for the energy ranges 



considered suggests that a = 0.5 might also correspond 
to the thermodynamic hmit for the interacting cases. 



A. Theoretical models 

We want to derive the entropy for a system of non in- 
teracting electrons with a constant density of states. It 
is possible to obtain analytically the entropy and other 
properties of this system for two models. Firstly, we 
study a continuous model using one-particle theory. Sec- 
ondly, we investigate a discrete model of equally spaced 
levels by counting directly the number of configurations. 

We first assume a continuous model with a constant 
density of states for the whole range of energies of inter- 
est, g{t) = go- The chemical potential fi and the Fermi 
level £/ are equal to zero by symmetry. At a given tem- 
perature T, the average occupation number of a state 
with energy e is: 



1 



1 -I- exp(e/T) ■ 



(5) 



So the mean energy (E) of the system, measured from 
the ground state energy, is 
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(6) 



where 9(e) is the step function, which takes into account 
the ground state energy, and V is the volume of the sys- 
tem. Eq. (P) allows us to change between temperature 
and mean energy. 

The entropy 5" of a non-interacting system can be cal- 
culated from the average site occupation through the ex- 
pression iQ 

S^J y5o|^.ln('^-lVln(l-n,)|de. (7) 

Integrating this expression and using Eq. (0) , one obtains 
the entropy as a function of energy: 



2TT^Vgo{E) 



(8) 



Eq. (||) is plotted as the solid line of Fig. 1. We can see 
that the numerical results do not match the theoretical 
one for the continuous case, Eq. m). To check the impor- 
tance of the discrete nature of our numerical systems we 
develop a discrete model. We now describe this discrete 
model for the number of configurations of non- interacting 
electrons with equally spaced one-particle energy levels. 
The dashed line in Fig. 1 represents log N{E) for this 
model. 

Each one-particle state is labeled according to its po- 
sition in the energy scale. The highest occupied level in 
the ground state is assigned to the zero of our scale. In 



the ground state, occupied levels correspond to negative 
numbers, while empty levels to positive numbers. A state 
k has an energy equal to kAe, where Ae is the constant 
energy spacing between levels, which we take equal to 1. 
The density of configurations n{E) can be decomposed 



as 



n(£;)=^n,(£;), 



(9) 



4=1 



where Ui (E) is the number of configurations with i exci- 
tations {i electrons above Fermi level and i holes below 
it) and with a total excitation energy equal to E. The 
excitation energy of a transition from level —m (with 
TO > 0), to level Z (^ > 0), is equal to / — to, recall that 
Ae = 1. For one electron jumps, ni{k) is simply 



ni{k) = k, 



(10) 



since we have k different possibles jumps of length k. 

For i — 2 (two excitations) the possibles jumps can be 
from £_„ to Ej and from e_p to e/, with the conditions 
"■ 7^ Pi J 7^ ^ (from two different states to two other 
different states), n,p > and j,l > 0. The total energy 
is the sum of individual energy 



k={j-n) + {l-p). 
We rewrite equation ( [ll| ) as 

k = j + I - {n+p). 



(11) 



(12) 



n2(fc) is, then, the number of possible combinations of 
j,l,n and p verifying expression (|l2|). Using the function 
Qi explained below and ( [12|) we arrive at 



fc-i 



n2ik) = }Q2U)Q2{k-j + 2). 



3=3 



(13) 



For an arbitrary i, the number ni(k) of i transitions 
with a total energy k becomes 



fe-i('i-l)/2 

ni{k)^ ^ Qz{j)Qz{k - j + i) 

j=i{i+l)/2 



(14) 



The integrated density of configurations N{k) is 

k 

N{k) = l + Y,n{i), (15) 

i=l 

where n{i) is given by Eq. (Q) and where we have added 
1 to count the ground state. 

We have to obtain the number Qi{m) of possible com- 
binations of i integers whose sum is equal to m, and 
subject to the constraint that they must all be different 
from each other and from zero. When i = 1 we have the 
trivial case 



Qi{m) = 1, 



(16) 



independently of to. It is easy to check that for i = 2 
that function becomes 



Q2{ni) = Int 



TO — 1 



(17) 



where Int (to) truncates to to its integer part. When 
i > 2, the function Qi{m) is defined in a recursive form 
as 



i{m) 



lnt{m/i) — l 



i{m-ij). 



(18) 



A hne connecting the points given by the integrated 
density of configurations N{E) of our discrete model, ob- 
tained from Eqs. (^, (|) and (|l|) (where E =k- Ae), 
has been plotted in Fig. 1 (dashed line) as a function of 
energy. We can see that our predictions with the discrete 
model fit very well our numerical results for systems with 
no interactions (solid dots). Furthermore, the number of 
configurations of the discrete model tends to the num- 
ber of configurations of the continuous model, Eq. (g), 
represented by the solid line in Fig. 1. 



same figure we also represent the results obtained by nu- 
merical calculations, for a non interacting system (solid 
dots) and for coulombic systems (triangles correspond to 
465 sites and empty squares to 899 sites) . The results for 
non-interacting systems agree fairly well with the discrete 
model and, using the normalized energy, are independent 
of size. On the other hand, the mean number of electrons 
above the Fermi level is larger for systems with interac- 
tions due to the electronic correlations, confirming our 
results of the previous section. It is interesting to note 
that the results for the interacting case are also indepen- 
dent of the size of the system when plotted, as in Fig. 2, 
using the normalized energy E, Eq. (Q) at least for the 
sizes considered. Both curves for Coulomb interacting 
systems fit to an expression of the form 



= C£;° 



(22) 



with an exponent a roughly equal to 0.25, half the value 
for non-interacting systems. 



10 



IV. DIFFERENCE IN OCCUPATION 

An interesting quantity is the number of electrons that 
have to jump from the ground state to reach a given 
configuration, i.e., half the Hamming distance between 
this configuration and the ground state. We consider the 
average {i) of this quantity over small energy intervals 
and disorder realizations. For a non-interacting system 
this magnitude corresponds to the average occupation of 
states above the Fermi level, which for a (continuous) 
constant density of states is given by: 



A 
V 
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1 + cxp(e/r) 



de = g^T\n2. 



(19) 



We use Eq. (pi) to rewrite {i) as a function of the mean 
energy, instead of temperature T, (£'), 
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FIG. 2. Number of electron above the Fermi level as a func- 
tion of energy. Solid dots correspond to a non-interacting 
system. Triangles correspond to a system of 465 sites and 
empty squares correspond to a system of 899 sites, both of 
them with Coulomb interactions. The dashed line is given by 
the discrete model and the solid line by the continuous one. 



We can also calculate this quantity using our discrete 
model, explained before, and we get 



(z)=iV-i(i^)(^5]jn,(z)|, 



(21) 



where N{E) is given by Eq. (|T^) and ni{E) by Eq. (14). 
In Fig. 2 we plot (i) given by Eq. (EOh (solid line) and 
Eq. (|2l]) (dashed line) as a function of energy. In the 



V. EFFECTIVE NUMBER OF ELECTRONS 

PARTICIPATING IN LOW-ENERGY 

TRANSITIONS 

Let us consider that our system is in any of the N{E) 
configurations up to energy E with equal probability and 
let us study how the site occupancy depends on the in- 
teractions present. For the non-interacting system all 
configuration are related to at least another lower energy 



configuration by a single electron jump. This is no longer 
true for interacting systems. We can generalize this idea 
through the concept of the Hamming distance between 
two configurations, which is defined as the number of sites 
with different occupation in both configurations. In the 
non-interacting case, for any given configuration there 
is always at least one lower energy configuration whose 
Hamming distance is equal to 2. We note that in the pre- 
vious section we calculated the Hamming distance to the 
ground state measured in the number of electron jumps 
needed to reach this state from a given configuration. 
Now we try to estimate for interacting systems the (av- 
erage) minimum Hamming distance from a given configu- 
ration to any lower energy configuration. Obviously this 
quantity will strongly affect the structure of phase space 
and in particular will influence the number of sites whose 
occupation changes in the N{E) configurations. 

For non-interacting systems following a Fermi-Dirac 
statistic, the importance of the number of sites whose 
occupation changes is reflected in the fact that, in the 
thermodynamic limit, the entropy S = \nN{E) can be 
expressed as a sum over the sites: 

n 

Sn = -J2 ((^*) 1^<*') + <1 - **) li^(l - ^*)) ' (23) 
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FIG. 3. Siv/S^^' as a function of ln(iV(S)), where N{E) is 
the number of low-energy configurations considered for sys- 
tems witli short range interactions of 899 sites (solids tri- 
angles) and systems with Coulomb interactions of 899 sites 
(empty dots), 465 sites (diamonds) and 248 sites (empty 
squares) . 



where (si) is the average occupation number at site i 
and n is the number of sites in the sample (note that 
sites with either {si) = or 1 do not contribute to Sn)- 
In general, for an interacting system, Sn does not cor- 
respond to the real entropy. A very simple interacting 
case where it is possible to obtain a similar relation is a 
system whose elementary excitations consist of simulta- 
neous jumps of q electrons. If these excitations are inde- 
pendent, the number of sites contributing to the sum in 
Eq. ( p3| ) is q times the number for the non-interacting 
case. Thus, Sn = q\nN{E) is proportional to the real 
entropy. We could propose Sn/ In N{E) as a rough quan- 
tity to estimate q but it is affected by the flnite size ef- 
fects observed in Fig. 1. To avoid this flnite size effects. 
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we propose the quantity Sn/S}^ , where Sj^ is Sn for 
a non-interacting system, as a measure of the number 
of electrons participating in a low-energy transition. In 
Fig. 3, we plot Sn/Sj^ as a function of \ii(N{E)) for a 
system with short-range interactions of 899 sites (solids 
triangles), and for three sizes of a coulombic system: 899 
sites (empty dots), 465 sites (diamonds) and 248 sites 
(empty squares) . Sn depends on the size of the system 
when interactions are present. For systems with short 
range interactions this dependence is much weaker and 
we only plot one size in Fig. 3 for the sake of clarity. 
These results nicely show the importance of correlations 
in the low-energy conflgurations of interacting systems. 



VI. VALLEYS IN CONFIGURATION SPACE 

To finalize, we study some aspects of the valley struc- 
ture of the space of low-energy configurations for systems 
with Coulomb interactions. In this space we can define 
a separation (which does not have the metric properties 
of a distance) between configurations / and J, with en- 
ergies Ej and Ej, respectively, through their transition 
rate ujij, taken to be equal to (l): 



w/j 



— exp I -2 2J rij I a \ exp <^ 



Ej — El 



(24) 



for Ej > Ej, and without the second exponential for 



the opposite case. In Eq. (24), tq is the inverse phonon 
frequency, of the order of 10~^^ s, a is the localization 
radius, and ^ r^ is the minimized sum of the hopping 
lengths of the electrons participating in the transition. 





FIG. 4. Relaxation-trees in configuration space. Each dot 
represents a configuration and each line the fastest transi- 
tion from the upper configuration, vertical distances repre- 
sent their energies, being their mutual horizontal distance the 
logarithm of the transition time. The bottom of the valley is 
a different configuration for each case represented. The con- 
figurations considered as bottom of the valley are the ground 
state (figure at the bottom), the second configuration and the 
eighteenth configuration (figure at the top). 

A way to study the valley structure is through the re- 
laxation processes. More specifically, we assume a very 
small temperature so that from a given excited config- 
uration we always go downwards in energy and we con- 
sider only the fastest transition rate at each step. Thus, 
we associate each configuration with a lower energy one 
and we construct in this way relaxation-trees. In Fig. 4 
we show three different relaxation-trees in configuration 
space for the same sample. All transition jumps consid- 
ered are plotted as lines where the horizontal distance is 
proportional to the logarithm of the transition time and 
the vertical one to the jump in energy. The bottom val- 
ley of the three cases represented are (in ascending order) 
the ground state, the second and the eighteenth configu- 
ration. At a qualitative level at least we can appreciate 
certain self-similarity between the trees. 

Each metastable configuration, whose fastest transi- 
tion downward in energy is carried out by the simulta- 
neous jump of two or more electrons, corresponds to the 
bottom of a valley. It can be characterized by the num- 
ber of electrons, i, participating in the fastest transition 
downward in energy. Its corresponding valley is the set 
of configurations that relax to it. For each characteristic 
number, i, let us define the integrated density of valleys, 
Ni{E), as the number of i-valleys up to an energy E. 
The standard integrated density of configurations is, ob- 



viously, No{E) = X]i=i ^i{E)j where in order to use a 
closed notation we denote by iVo to the standard inte- 
grated density of configurations and by Ni to the inte- 
grated density of configurations whose fastest transition 
is a one electron jump. We have numerically checked that 
these integrated densities fit fairly well an expression of 
the form 



N,{E)^Ce^Y>{{E/Ear}, 



(25) 



where C and Eo depend on the characteristic number of 
electrons. On the other hand, the exponent a does not 
seem to depend on the characteristic number of electrons 
although it presents very large fluctuations from sample 
to sample. For the range of energies studied this expo- 
nent is roughly 0.7, the exponent for the integrated den- 
sity of configurations Nq{E) at the energy interval con- 
sidered. We know from a previous section that a = 0.5 in 
the thermodynamic limit. The proportion of metastable 
states with respect to normal configurations decreases 
with increasing energy. In Fig. ^, the density of configu- 
rations Nq{E) (empty squares) and the density of valleys 
stables against 2 (solid dots), 3 (solid diamonds) and 4 
(solid triangles) simultaneous transitions are plotted as 
a function of energy, E. All curves fit fairly well the ex- 
pression N{E) = Cexp{E/Eo)", with different C and Eq 
but with an exponent a very close to 0.7 in all cases. 
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FIG. 5. Integrated densities of configurations Ni{E) versus 
energy E. Empty squares corresponds to the integrated den- 
sity of configurations, i = 0, while the other curves correspond 
to integrated densities of metastable configurations with dif- 
ferent characteristic number of electrons: 2 (solid dots), 3 
(solid diamonds) and 4 (empty triangles). All curves fit fairly 
well the expression N{E) = C exp{E / Eo)" , with different C 
and Eq but with an exponent a close to 0.7 in all cases. 



VII. CONCLUSIONS ping and Related Phenomena, edited by O. Millo and Z. 

Ovadyahu, p. 35, (Racah Institute, Jerusalem 1995). 
We have calculated the integrated density of config- [12] K. Huang, in StattsUcal Meehamcs Wiley, (1963) 
urations for systems with no interactions, short range 
interactions and Coulomb interactions. A coniparation 
of our numerical results with a continuous theoretical 
model shows important finite size effects. However, our 
data overlap fairly well with a discrete model in the non- 
interacting case. The finite size effects found in the non- 
interacting case could serve to estimate these effects in 
interacting systems. We have more deeply studied the 
structure of phase space analyzing other quantities like 
(?) (half the average Hamming distance from a given con- 
figuration to the ground state), Sm (the "site occupation 
entropy" see Eq. (p3)) and the valleys in configuration 
space. For a given energy, (i) strongly depends on the 
type of interactions while, for the systems considered, it 
does not depend on the system size. The "site occupation 
entropy" , Sm , increases with the strength of interactions 
and, for a long range one, with the size of the system. 
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